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Abstract 

Bistable structures associated with non-linear deformation behavior, exemplified by the Venus 
flytrap and slap bracelet, can switch between different functional shapes upon actuation. Despite 
numerous efforts in modeling such large deformation behavior of shells, the roles of mechanical and 
nonlinear geometric effects on bistability remain elusive. We demonstrate, through both theoretical 
analysis and table-top experiments, that two dimensionless parameters control bistability. Our 
work classifies the conditions for bistability, and extends the large deformation theory of plates 
and shells. 

PACS numbers: 46.25.-y, 02.40.-k, 62.20.F-, 02.30.Oz 
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Mechanical structures exhibiting multiple s table morphologies arise in a variety of physi- 
cal systems, both natural [1] and synthetic j2"5|. For example, the leaves of the Venus flytrap 
can be triggered to snap-through from an open to a closed state in a fraction of a second 
to capture insects Indeed, multi-stable components have promising applications in mi- 
cropumps, valves, shape-changing mirrors in adaptive optical systems, mechanical memory 
cells [5|], artificial muscles, bio-inspired robots Ig], deployable aerospace components [31], and 
energy harvesting devices 



Multi-stability of plates and she 



many studies^ over the years 
structures 



Is associated with non-linear deformation has inspired 
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16| . while related bifurcation phenomena in slender 



17l |. such as DNA [18|, ll9| and plant tendrils [20|, have also been long investigated 



2l| . Existing models of shell structures utilize geometric assumptions or scaling arguments 



that greatly simplify analysis, leading to solutions that have significantly improved our un- 
derstanding of bistability. Nonetheless, these geometric assumptions are typically only valid 
for small angle deformations and often do not strictly satisfy the mathematical conditions 
for geometric compatibility ISj, especially in large deformation cases. Scaling theory has 
been employed by Forterre et al. ij, who proposed a minimal model to explain some of the 
bistable features of Venus flytraps, with a dimensionless deformation energy postulated in 
lieu of a "first-principles" derivation, and Armon et al. [3], who investigated the competi- 
tion between bending and stretching deformation within incompatible thin sheets. In the 
latter work, a dimensionless width was proposed to be the key bifurcation parameter, with 
the threshold value obtained through experiments. While this study is elegant, it did not 
yield quantitative predictions for the bifurcation threshold; furthermore, the presence of an a 
priori unknown quantity (curvature) in the bifurcation parameter further restricts its utility. 
Therefore, for large deformations of shells, a more general elasticity theory framework, that 
explicitly includes the non-uniform bending curvature and mid-plane stretching, is required. 

In this work, we propose an analytically tractable, reduced-parameter model that in- 
corporates geometric nonlinearity and accounts for the competition between bending and 
in-plane stretching involved when a planar sheet is mapped into a non-developable shape 
via misfit strains. We employ this model to perform a comprehensive examination of bista- 
bility, leading to unique predictions for bifurcation thresholds and bistable morphologies 
that are quantitatively validated with a series of table-top experiments. More generally, our 
approach provides a quantitative means to investigate the role of geometric nonlinearities 
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FIG. 1: (A) Schematics: two thin latex rubber sheets (blue and yellow) were pre-stretched along 
perpendicular directions and bonded to a much thicker elastic strip. When released, the bonded 
multilayer sheet deforms into one of the following shapes. Here, H denotes the overall thickness, 
while angle (p measures the misorientation of the principal axes of curvature (ri and from the 
principal geometric ones (d^; and dy, respectively). (B) A saddle shape for a small, thin square 
sheet {H = 1.5mm, L = W = 24.0mm). (C) A saddle shape for a thick square sheet {H = 2.5mm, 
L = W = 48.0mm). (D) A stable, nearly cylindrical shape (curving downwards) for a thin, wide 
strip (H = 1.5mm, L = W = 48.0mm). (E) The other stable, nearly cylindrical configuration 
(bending upwards) for the same sheet as in (D). The nearly cylindrical shape smoothly transitions 
into a doubly curved shape near the edges. 



in a broader setting, including morphogenesis 



22h25| . film-on-substrate electronics 



26-28], 

and spontaneous bending, twisting and buckling of thin films and helical ribbons [3014331. |35|. 

We employ a simple experimental demonstration of the morphological transition from 
mono- to bi-stable states when the dimensions change (either the width increases or thickness 
decreases), with all other parameters fixed. To this end, two pieces of thin latex rubber 
sheets (thickness Hi, length L, and width W < L) are pre-stretched by an equal amount 
and bonded to an elastic strip of thicker, pressure-sensitive adhesive (thickness H2) [36^ 
as shown schematically in Fig. [1] A, such that the total thickness of the bonded strip is 



H = 2Hi + H2. The pre-strains along the two lateral directions were chosen to be equal to 



3 



0.28 in all of the experimental results presented. When released, the initially flat, bonded 
laminate deforms into either a saddle shape (Fig. [T]B and C) or one of two nearly cylindrical 
configurations (Fig. [T]D and E), driven by misfit strains. More specifically, when the strip 
is sufficiently narrow or thick, the equilibrium saddle shape is unique; while if the strip is 
wide and/or thin, it will bifurcate into one of the two nearly cylindrical shapes. 

As the starting point of our theoretical approach, we assume that an originally fiat 
elastic sheet of H W < L, with a rectangular cross-section and principal geometric 
axes along its length (d^), width (d^^), and thickness (d^), deforms into part of a torus 
shape with variational parameters ni and k,2 {i^i = ^/Ri, i = 1;2) along principal direc- 
tions ri and r2 (which coincide with the geometric axes) as shown in Fig. 2. A torus can 
describe a broad range of possible morphologies, such as a saddle, cylinder, and sphere; 
selected by tuning two geometric shape parameters. Instead of assuming small deforma- 
tion as in classical lamination theory, we allow for large deformations with the associ- 
ated geometric nonlinearity in the small elastic strain limit (in the thick sheet); that is, 
we only require that HiW <^ 1. Finally, mapping the sheet onto the surface of a torus 



(or part thereof) facilitates the explicit construction of the strain tensor components |39| : 
7^^(t, z) = [cos{K2t)z]/ {1/ki + [cos(K2t) - 1]/k2} + [cos(/t2t) - 1]/«i/k2, ^yy{t, z) = K2Z, and 
■yzz(t,z) = —uijxx + lyy)/i)- ~ ^) 1 where t and z denote the arclength along Y2 and the 
distance from the mid-plane of the shell. 

The strip is subjected to effective surface stresses (due to the thin, misfitting top and 
bottom layers) of the form f+ = ~f2^y ® dy and f ^ = fid^ ® d^^ along the top and bottom 
surfaces ("(g)" denotes the tensor direct product), respectively. This is equivalent to the 
case in which a strip is subjected to only one surface stress f ~ = /id^ (g) d^ + fod,, ® d^^ 
on the bottom surface, where only the bending mode of deformation is of interest 321 • -^^^ 
conciseness, we assume that the principal axes of bending, ri and r2, coincide with the 
principal axes of the surface stresses, i.e., the geometric axes, d^,. and d^ (justified below). 

Continuum elasticity theory gives the potential energy density per unit length of the strip 



as n = [f • l\z=~H/2 + j"H,2 ^1 ■■ C : jdz 

stiffness tensor. By employing the expressions for the strain components and expanding in 
kW <^ 1 (and noting that H <^ W), it is straightforward to show that 

U = Us + Ub + Ug + 0{EH'^W^k\ EHW'^k^), (1) 



H/2 1 



dt, where C denotes the fourth-order 
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FIG. 2: Deformation of an originally flat elastic strip into a section of a torus with variational 
parameters Kj (kj = 1/Ri, i = 1,2). Here, ri and r2 denote the principal bending axes and t 
denotes the arclength measured from the circle of radius Ri. The red dotted line at t is shorter 
than the red solid line, indicating an extra strain due to geometric constraints. 

where lis = "(/i^^i + f2i^2)WH/2 denotes the contribution from the surface stresses, Ilf, = 
{kI + Kl + 2iyKiK2)EW /2A{l-u^) is the bending energy, and Ug = EHW^{KiK2f /Q40{1- 
z/^) is the stretching energy arising from the geometric nonUnearity associated with non-zero 
Gauss curvature. In general, terms associated with bending deformation are of the form 
~ EW K,"^ {nWy^ , where n = 0, 1, 2, 3, while those associated with Gaussian curvature 
are of the form ~ EHW^k,'^{k,W)'^'^, where m = 0,1,2,3,.... Also, we note that in the 



special case of spherical bending (i.e., ki = K2 = 1/p), 11^, = EHW^J 



40 



640(1 - z/2)p^] to 



4l|. 



leading order, in agreement with classical Foppl-von Karman theory 

Since Ilf, ~ EWH^k,^ and Ug ~ EHW^k,'^, a simple scaling analysis suggests that when 



W <^ ^jHjK^ life ^ Ilg, and the nonlinear geometric effect becomes negligible. That is. 



bistability is controlled by a dimensionless geometric parameter 77 = W \/ k/ H [13], where 
the curvature k = max{|/ti|, \ n2\}- In the limit rj 1, Ug can be ignored, and the laminate 
develops into a saddle shape, as shown in Figs. [T]B and C (when /i = — /2). In this case. 
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applying the stationarity conditions [dH/dKi = 0, i = 1,2) to Eq. ([T]) and neglecting Eg, we 



recover the analytical results in [32|: ni = 6(/i — z//2)/-Eif^ and K2 = 6(/2 — i'fi)/EH^, 
with n* = —3{fi — 2z//i/2 + f2)/2EH. On the other hand, in the limit where 77 — )■ 00, 
the system bifurcates into two stable, nearly cylindrical configurations, for which one of the 
principal curvatures becomes very small over the interior of the shell, as shown in Figs. [T] D 
and E and the supplementary movie. 

The parameter t] is essentially the same as the one proposed by Armon et al. [iS^ . While 
rj quantitatively captures the effect of the geometric nonlinearity, it contains an unknown 
parameter, the curvature - significantly restricts its utility. This deficiency is remedied 
as follows. Since k ~ f /EH^ (where / = max{|/i|, I/2I}), we define a new dimensionless 



parameter, C = f /EHW/H , which involves only a priori known material and surface 
stress parameters. We demonstrate that this new parameter arises naturally by considering 



the case /i = — /2 relevant for our experiments and those of Armon et al. 13|, and derive 
an analytical expression for the bifurcation threshold. 

Applying the stationarity conditions to Eq. ([T]) and setting /i = — /2 yields (ki + 
K2) [/^i/^2 + 80(1 + i')H'^ /SW^] = 0. Therefore, either ki + K2 = which corresponds to 
the saddle-like shape, or hik,2 = —80(1 + u)!!'^ /SW^, which, when real solutions exist, 
yields the bifurcated solutions. A bifurcation occurs when both ki + K2 = and = 



-80(1 + iy)H^/3W^ are satisfied, implying that ^ = -K2 = k = -^80(1 + iy)/3H/W^ 
W^n/H = r/c = [80(1 + z/)/3]^/^ Importantly, our prediction for the bifurcation 
threshold {rjc = 2.51 for u = 0.49) is in excellent agreement with the experimental result, 



rjc ~ 2.5, found by Armon et al. 13|. Moreover, stationarity conditions at the bifurcation 
imply that ki-K2 = 2k = 6f*{l-u^)/EH^, or k = 3f*{l-iy^)/EH^. Thus, the bifurcation 
threshold can also be expressed as 

.1/4 



80 



27(1 - z/)(l - z/2 



(2) 



EH H ^3(1-1/2) 

Equation (|2]), which is the central result of the analysis, quantitatively expresses the funda- 
mental condition for bistability in terms of a dimensionless parameter which only depends 
on the driving force (surface stress), material properties, and ribbon dimensions. Although 
the bifurcation formally occurs aX ( = (c, since the shape change only takes place gradually 
as ( increases, we expect the bistable behavior to manifest only when ( ^ (c- 

In comparing theoretical predictions with experiments, we find that in both Figs. [TJ3 








FIG. 3: Bistability map in the (/i,/2) space. In (A), the red shadmg indicates bi-stable regions 
while no shading indicates mono-stability. The red arrows denote transitions from mono-stable to 
bistable states. The total strain energy versus the orientation of the misfit axis for (B) /2//1 = \/3, 
(C) /2//1 = 1, (D) /2//1 = 1/^/3, (E) /2//1 = -1/^/3, (F) /2//1 = -1, and (G) /2//1 = -^/3. 

and [IP, ( is comparable to Cc ~ 1-5 (in B C = 1.5; in C C = 1-4), and the sheets adopt 
saddle shapes, as expected. On the other hand, increasing ^ to 3.1 leads to bistable behavior 
characterized by cylindrical shapes, as shown in Figs. [Tp and [TJ3, respectively. Moreover, 
for the cylindrical shapes, the measured radius of curvature (13.1 ± 0.5mm) is in reasonable 
agreement with the theoretical prediction EH'^/[Qf{l — z/^)] ^ 13.8mm. 

We now relax the condition fi = —j^ and consider the role of surface stresses on me- 
chanical instability in more general terms. In the C — > 00 (77 — > 00) regime, the geometric 
nonlinearity requires that either /ti — )■ or ^2 — ?■ away from edges. Without loss of gen- 
erality, we assume that Ky = ^ 0. In order to examine stability, we first note that the 
principal bending axes ri and r2 do not necessarily coincide with the geometric axes and 
dj^, but instead form an angle as shown in Fig. [TJ\. In this case, the potential energy per 
unit length of the strip becomes H = -W H^hC^ ^ JiS'^)^^!'!^ EE'^W k\I2^{\- v^), where 
C = cos (f) and S = sin 0. Minimizing 11 with respect to both Hi and (p yields two sets of so- 
lutions: (1) 0* = with Ki = 6(1 - z/2)/i/Eif2 or (2) (p* = tt/2 with ki = 6(1 - u^)f2/EH^. 
Whether an extremum configuration (0 = or = 7r/2) is locally stable depends on whether 
n"(0*) = d^U/d(f)\=^. > 0. To this end, we write n"(0) = -6(1 - u^){f2 - fi)fiW/EH oc 
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- 1) and n"(7r/2) = -6(1 - iy')ifi - f2)f2W/EH oc - /3), where /5 ^ h/fi 

denotes another key dimensionless parameter. Clearly, whenever /3 < 0, both n"(0) > and 
n"(7r/2) > 0, implying the emergence of two mechanically stable configurations with one 
configuration in general being energetically more favorable than the other. Interestingly, 
for /3 = — 1, the two stable configurations are degenerate and have the same energy in the 
absence of edge effects. On the other hand, when < /3 < 1, n"(0) > while n"(7r/2) < 0, 
implying that the 0* = configuration is stable, while the (p* = 7r/2 one is mechanically 
unstable. For /3 > 1, the stabilities of the two extrema are reversed, with the 0* = 7r/2 
configuration now displaying mechanical stability. Finally, for th e sp ecial case /3 = 1, n*(0) 
is constant, implying that the shell is in a neutrally-stable state 42|. In this case, the lami- 
nate can, in principle, curl into a nearly cylindrical shape with arbitrary orientation (at least 
in the absence of edge effects), and can also easily transition between these shapes drawn 
from a continuous family of neutrally-stable configurations. The different possibilities are 
illustrated in Fig. |3l 

Finally, when edge effects are taken into account, even when /3 = ±1, an asymmetry in 
dimensions {L ^ W) can result in the breaking of the degeneracy between the minimum 
energy configurations. As shown in Fig. [1] E, the cylindrical shape in the ( ^ (c limit 
smoothly transitions into a more saddle-shaped one with finite Gauss curvature over a 



length scale ~ Wc ~ QH^yEH/ f as the ends of the cylinder are approached. Within 
this transition region, the total potential energy per unit area can be approximated by 
HI = — 3(/f — 2z//i/2 + /|)/2£'if, while for the same area in the cylindrical part, the energy 
per unit area is = -3(1 - i''^)fl/2EH. Since Yi\ < Yi^ for /3 = ±1, the edges reduce the 
overall deformation energy, and the most effective reduction is obtained for maximizing the 
overall edge length. Thus, the laminate will curve along the long direction, in agreement 
with the recent results of Alben et al. j43|. 

In summary, large shell / plate deformation with geometric nonlinearity is treated through 
a novel, analytically tractable theoretical framework which combines continuum elastic- 
ity, differential geometry and stationarity principles. Two key dimensionless parameters 
(C = \J f / EHW/H and P = /2//1) were shown to govern structural bistability. It is note- 
worthy that C only involves the driving force (surface stress), material properties, and ribbon 
dimensions, in contrast to the geometric bifurcation parameter rj = W^/ k/H introduced by 



Armon et al. 



13|, which involves an a priori unknown parameter (curvature). On the one 
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hand, whether a structure with hnear elastic properties can exhibit bistabihty depends on 
its dimensions and curvatures (C ^ (c, or rj ^ rjc) due to nonhnear geometric effects. On the 
other, even when the necessary geometric condition is satisfied, the existence of bistabihty 
still depends on a mechanical parameter (/3 < 0). Our theoretical analysis also predicts the 
lifting of ground state degeneracy (when (3 = ±1) due to edge effects [43|. The non-linear 
geometric effects on bistabihty have been verified by table-top experiments, and our theo- 
retical predictions for the bifurcation threshold are also in agreement with the experimental 
results of Armon et al. [13| . 

In a broader sense, our approach provides a means to quantify and understand the role 
of nonlinear geometric effects in multi-stable structures in a wide range of two-dimensional 
natural and engineered systems, such as morphogenesis and mechanical instability in biolog- 
ical systems, film-on-substrate electronics, and spontaneous coiling and buckling of strained 
thin films and ribbons. Via reverse-engineering, it will also facilitate the design of multi- 
stable functional structures, from bio-inspired robots with smart actuation mechanisms to 
deployable, morphing structures in aerospace applications. 
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